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ABSTRACT 

We present new near- and mid-IR observations of 19 Class I/flat-spectrum young 
stellar objects in the nearby p Ophiuchi {d = 125 pc) and Serpens {d = 310 pc) 
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dark clouds. These observations are part of a larger systematic infrared multiplicity 
survey of Class I/flat-spectrum objects in the nearest dark clouds. We find 7/19 
(37% lb 14%) of the sources surveyed to be multiple systems over a separation range 
of ~ 150 - 1800 AU. This is consistent with the fraction of multiple systems found 
among older pre-main-sequence stars in each of the Taurus, p Ophiuchi, Chamaeleon, 
Lupus, and Corona Australis star-forming regions over a similar separation range. 
However, solar-type main-sequence stars in the solar neighborhood have a fraction 
approximately one-third that of our Class I/flat-spcctrum sample (11% it 3%). This 
may be attributed to evolutionary effects or environmental differences. An examination 
of the spectral energy distributions (SEDs) of the SVS 20 and WL 1 binaries reveals 
that the individual components of each source exhibit the same SED classifications, 
similar to what one typically finds for binary TTS systems, where the companion of a 
classical TTS also tends to be of the same SED type. 

Subject headings: binaries: close — stars: formation — stars: pre-main-sequence 

1. Introduction 

The past two decades have witnessed substantial progress in the study of star formation and 
its aftermath, the birth of planetary systems. Observations have led the way, with much of the 
credit due to both spaccborne missions, such as the Infrared Astronomical Satellite (IRAS), the 
Infrared Space Observatory (ISO), and the Hubble Space Telescope (HST), and to ground-based 
infrared (IR) imaging surveys from NOAO/SQIID, the NASA Infrared Telescope Facility (IRTF) 
and other facilities. Advances in detection and instrumentation have been matched by equally 
exciting theoretical developments. Researchers have delineated key features of a young star's 
structure and the manner in which it affects its surroundings, both thermally and mechanically. 

As we broaden our focus from individual stars, we find that there is a universal tendency for 
these objects to form in groups, rather than as isolated entities. On the largest scales, IR surveys 
have revealed stellar aggregates still embedded in their parent molecular clouds (e.g., Zinnecker, 
McCaughrean, &; Wilking 1993). These must eventually become optically visible associations or 
bound clusters. However, the smallest entities created within such groups are generally not single 
stars. We have known for many years that the majority of solar-to-late-type field stars are binaries 
or multiples (Abt k, Levy 1976; Duquennoy h Mayor 1991; Fischer &; Marcy 1992), but only 
in the past decade have significant numbers of pre-main-sequence (PMS) T Tauri stars (TTSs) 
been surveyed for multiplicity. These recent IR surveys have shown that young low-mass stars 
have binary fractions that are greater than or equal to that of the field (e.g., Ghez, Neugebauer, 
k. Matthews 1993; Mathieu 1994; Simon et al. 1995). These pairs appear to be coeval, that is, 
at the same evolutionary age (Hartigan, Strom, k, Strom 1994; Brandner & Zinnecker 1997). 
Recent observations at millimeter continuum wavelengths also find multiple sources at the origins 
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of extended molecular outflows and optical jets (Looney, Mundy, &; Welch 1997). The current 
findings set this fact in its proper evolutionary context. Most stars form as part of a binary 
system, rather than find their partners later in life. 

The large binary frequency in visible, PMS stars prompts us to investigate even younger 
systems. Recent sub-arcsecond surveys of young stars have almost exclusively focused on PMS 
TTSs and have virtually ignored younger systems with Class I/flat-spectrum energy distributions. 
Thus we know very little about the multiplicity of self-embedded young stars. Is multiplicity the 
rule even among more deeply embedded objects, those still forming out of their parent clouds? 
What is their binary frequency, and how does this compare to that of PMS stars? What are the 
energy distributions of their components? New multiplicity surveys of known embedded low-mass 
stars must be conducted to address these questions. 

We have initiated the first systematic near-to-mid-IR multiplicity survey of ~ 100 Class 
I/fiat-spectrum young stellar objects in six nearby (d <^ 300 pc) dark clouds. Our survey has 
several purposes: 1) to study the binary fractions and separations of protostars; 2) to diagnose 
the luminosities and evolutionary states of their components; and 3) use our near-IR photometry 
(in conjunction with existing IRTF and Keck spectra) to place protostars in Hertzsprung- Russell 
(H-R) diagrams to diagnose their stellar properties. 

All of the sources in our survey were selected such that they have either Class I (protostellar) 
or fiat-spectrum (slightly more evolved than Class I) spectral energy distributions (SEDs; see 
Adams, Lada, & Shu 1987; hereafter ALS87, 1988; Lada 1987). These objects are the direct 
evolutionary progenitors of the well-studied TTSs, and many are reasonably bright in the near-IR 
[vciK ~ 10). Although these objects have peak fluxes at far-IR wavelengths, there are no far-IR 
high-resolution observatories (0.85 - 3.0 m for SIRTF - SOFIA) currently operating. Also, 
millimeter wavelength observations of the very few known Class I binaries show that their envelopes 
can be so large that they overlap, making it difficult to discern the individual components at these 
long wavelengths (Looney, Mundy, k, Welch 2000). 

In this contribution, we present the initial results of our survey. Specifically, we discuss our 
near- and mid-IR observations of 19 Class I/flat-spectrum objects in the p Ophiuchi [d = 125 pc; 
Knude & Hog 1998) and Serpens (d = 310 pc; de Lara, Chavarria-K., &: Lopez- Molina 1991) dark 
clouds. We discuss the observations and data reduction procedures in §2. In §3, we present the 
results of our survey, and discuss the results in §4. We summarize our primary results in §5. 

2. Observations and Data Reduction 

Two different sets of infrared observations were obtained for this survey. The near-IR 
Ji^i^L-band (1.25, 1.65, 2.2, and 3.5 /xm) observations of each source were obtained with the 
NSFCAM 256x256 InSb facility array camera (Rayner et al. 1993; Shure et al. 1994) on the 
NASA IRTF 3 m telescope on Mauna Kea, Hawaii. These data were supplemented with mid-IR 
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A?'-band (10.8 /um) observations taken with MIRLIN, JPL's 128x128 pixel Si:As camera (Ressler 
et al. 1994), on both the Palomar Observatory 5 m and Keck II 10 m telescopes. A summary of 
all observations is shown in Table 1. 



2.1. NSFCAM-IRTF Near-IR Observations 

The J, H, K and L-band observations of each source were made during the period 2001 July 
08 - 10. For our study, we used a plate scale of 0'.'148 pixel"^ with a corresponding field of view of 
approximately 38"x38". The average FWHM for all observations was ~ 1". 

Each source was observed in a five point dither pattern (a 2 x 2 square with a point at 
the center) with 12" offsets between the corners of the square. At each dither position, the 
telescope was nodded to separate sky positions 30" north of each target observation. Typical total 
integration times ranged from 10 seconds to 15 minutes for the JHK data, and 60 seconds at 
L-band. 

All JHKL data were reduced using the Image Reduction and Analysis Facility (IRAF)^. 
The individual sky frames were normalized to produce flat fields for each target frame. All target 
frames were processed by subtracting the appropriate sky frames and dividing by the fiat fields. 
Finally all target frames were registered and combined to produce the final images of each object. 

2.2. MIRLIN-Keck-Palomar Mid-IR Observations 

All mid-IR A^'-band (Aq = 10.78 /nm, AA = 5.7 fim) data were obtained as part of a larger, 
mid-IR imaging survey of the p Oph cloud core's embedded population (Barsony, Ressler, &; 
Marsh 2002). The plate scale and field of view of MIRLIN at Keck II are ^.'138 pixer^ and 17'.'7 
X 17'.'7, respectively. Corresponding values at the Palomar 5-meter are 0'.'15 pixel"^ and 19^'2 x 
19'.'2. For reference, the full-width at half-maximum of a diffraction-limited image at A'^-band is ~ 
0'.'25 at Keck II and 0'.'47 at the Palomar 5-meter. 

Data were acquired with traditional mid-IR chopping and nodding techniques. Specifically, 
the telescope's secondary mirror was chopped 8" in a north-south direction, at a rate of a few 
Hz; then the entire telescope was nodded 8" east-west in order to remove residual differences in 
the background level. Total on-source integration times were typically 24 - 25 seconds at each 
telescope for the program sources. On-source integration times consisted of several hundred to 
a thousand coadded chop pairs, with 5-6 msec integration times per frame. The raw images 
were background-subtracted, shifted, and co-added with our in-house IDL routine "MAC" 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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(match-and-combine) . 

2.3. Near-IR Source Photometry and Calibration 

Aperture photometry was performed using the PHOT routine within IRAF. An aperture of 4 

pixels in radius was used for all target photometry, and a 10 pixel radius was used for the standard 
star photometry. Sky values around each source were determined from the mode of intensities in 
an annulus with inner and outer radii of 10 and 20 pixels respectively. Our choice of aperture size 
for our target photometry insured that the individual source fluxes were not contaminated by the 
flux from companion stars, however they are not large enough to include all the flux from a given 
source. In order to account for this missing flux, aperture corrections were determined using the 
MKAPFILE routine within IRAF. The instrumental magnitudes for all sources were corrected to 
account for the missing flux. 

Photometric calibration was accomplished using the list of standard stars of Elias 
et al. (1982). The standards were observed on the same nights and through the 
same range of airmasses as the target sources. Zero points and extinction coefficients 
were established for each night. All magnitudes and colors were transformed to the 
CIT system using Mauna Kea to NSFCAM and NSFCAM to CIT photometric color 
transformation equations from http://irtf.ifa.hawaii.edu/Facility/nsfcam/mkfilters.html, 
http://irtf.ifa.hawaii.edu/Facility/nsfcam/hist/color.html and the NSFCAM User's Guide. The 
photometric accuracy for all observations is typically good to within it 0.10, 0.03, 0.03, and 0.06 
magnitudes at J, if, K, and L-band respectively. GY 224 has a J-band magnitude uncertainty of 
dr 0.32 mag, and the iJ-band magnitude uncertainty for EC 53 is dr 0.08. 

2.4. Mid-IR Source Photometry and Calibration 

Flux standards were a Corona Borealis (iV = +2.19) and 7 Aquilae {N = -0.78) for the Keck 
II run, and 7 Aquilae and Ophiuchi {N = +2.20) for the Palomar data. These were also used 
for airmass monitoring during each night of observing. 

Photometry for the standard stars was performed in 14-pixel radius apertures (corresponding 
to a 1''9 radius software aperture at Keck II and a 2'!'! radius aperture at Palomar). A straight-line 
fit to the instrumental magnitudes as a function of airmass for each night resulted in the 
determination of the zero-point offsets and airmass corrections. No airmass correction was 
required for the Keck II data, whereas typical airmass corrections at Palomar were of order 0.1 
- 0.2 mags/airmass. The photometric consistency between all standards during a given night's 
observing was typically of order ± 0.04 magnitudes. By adding the errors in the zero-point offsets, 
the airmass corrections, the aperture corrections, and the uncertainties in the magnitudes of the 
standards in quadrature and taking the square root, we estimate the total photometric accuracy of 
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the Keck data to be good to it 0.06 magnitudes at AT-band. The Palomar data have a photometric 

accuracy good to it 0.08 magnitudes at A'^-band. To convert these errors to Janskys, note that 
0.00 magnitudes at iV-band with MIRLIN corresponds to 33.4 Jy. Photometry for program 
sources was performed in a 5-pixel radius aperture (corresponding to 0'.'75 and 0'.'69 radius software 
apertures at the Hale 5-meter and at the Keck II 10-meter telescopes, respectively). Aperture 
corrections were derived from the flux standards for each night, and applied to each target object's 
instrumental magnitude before application of the zero-point calibration and airmass correction. 
In the case of close companions, the combined system flux was first determined from a software 
aperture chosen to be large enough to contain both objects. Subsequently, the relative photometry 
of each component was determined by fitting a known point-source calibrator (generally one 
of the flux calibrators observed that night) to the individual source peaks. The total flux was 
then divided amongst the components in the ratio determined by the relative point-source fltting 
photometry. 



3. Analysis and Results 

Five of 13 (38% ± 17%) of the p Oph sources (sepn. range ~ I"- 10'.'5 or ~ 125 - 1500 AU) 
and 3/6 (50% it 29%) of the Serpens sources (sepn. range ~ 1'.'5 - 5"or ~ 470 - 1560 AU) appear 
to be binary or multiple in nature. Thus, among the sources surveyed, we find 8/19 (42% it 15%) 
stars in the separation range 1"- 10'.'5 (~ 125 - 1560 AU) to be members of binary or higher order 
systems. The quoted uncertainties represent the statistical \/iV errors in the fraction of sources 
found to be binary or multiple. In Figures 1 - 8, we present our near and mid-IR images of each 
multiple source. IRS 43 was not detected at J-band, and we do not as yet have 10 ^m data 
for L1689 SN02 and EC 82/EC 86. Because the spacing of our on-chip chop/nod pattern was 
close to the component separations for GY 23/GY 21 and EC 92/EC 95, wc do not present 10 
micron images of these systems. The photometry was, nevertheless, recovered from the raw data. 
Contour levels for each image are listed in the respective flgure captions. EC 82 forms a possible 
wide binary with EC 86. However, given their separation of ~ 2700 AU at our adopted distance 
to Serpens of 310 pc, these objects may not compose a true bound binary system. We therefore 
exclude this system from our calculated CSF. Each source is individually discussed in detail in the 
next section. 

In Tables 2 and 3, we present the JHKL magnitudes, 10 /im fluxes (in Janskys), near-IR 
colors, and spectral indices for all surveyed sources. Spectral indices were calculated from a least 
squares fit to the 2.2 /im to 10 iim data. We do not as yet have MIRLIN A^-band fluxes for IRS 
63, IRS 67, VSSG 17, L1689 SN02, ISO 159, EC 53, and EC 82/EC 86. However, mid-IR ISO 
fluxes at 6.7 /xm and 14.3 /xm have been published for the p Oph single sources IRS 67 and VSSG 
17 by Bontemps et al. 2001. The spectral indices reported in Table 2 for these two sources were 
calculated using the ISO mid-IR fluxes. If neither MIRLIN nor ISO fluxes were available for a 
given source, no spectral index was calculated. 
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Seven sources in Tables 2 and 3 were not detected at J-band. This was not likely due to faint 
sources being smeared out in these relatively long (90 second) exposures, since this effect was 
not observed in other equally long J-band images in which faint sources were detected. Upper 
limits for the J-band magnitudes for these sources were determined by adding artificial stars 
to the respective J-band images, and counting the number of sources recovered by DAOFIND. 
Artificial stars were added at random positions to each image in twenty separate half magnitude 
bins with each bin containing one hundred stars. The twenty bins covered a magnitude range 
from 15.0 to 25.0. The artificial stars were examined to ensure that they had a similar FWHM of 
the point-spread function as the sources detected in other J-band images. Aperture photometry 
was performed on all sources to confirm that the assigned magnitudes of the added sources agreed 
with those returned by PHOT. All photometry agreed to within 0.1 magnitudes. DAOFIND and 
PHOT were then run and the number of identified artificial sources within each half magnitude 
bin was tallied. This process was repeated 20 times. Our final J-band magnitude limits, good to 
within lb 0.25 magnitudes, are listed in Tables 2 and 3. 

In Table 4, we list the separation (in arcsec) and position angle (measured with respect to the 
brightest source at i^-band; the southernmost source in each case) for the sources in our survey 
which were found to be multiple. Figures 9 and 10 show the SEDs for WL 1 and SVS 20. In each 
case, the SED of the northernmost component has been shifted down by 1.0 in the vertical axis 
for clarity. The p Oph source L1689 SN02 and the Serpens sources EC 82/EC 86 were also found 
to be binary, however no resolved 10 fim. fluxes are currently available for these objects, therefore 
SEDs were not constructed. Wc can obtain estimates for the SED classes of these objects by 
examining their location in JHK and JHKL color-color diagrams. In this respect, EC 86 is found 
to have near-IR colors indicative of a reddened background or Class III source (PMS star with 
no circumstellar material; ALS87). By contrast, the colors of both components of L1689 SN02 
and EC 82 show infrared excesses suggestive of Class II (TTS surrounded by a circumstellar disk; 
ALS87) objects or earlier. 

4. Discussion 

4.1. Notes on the Multiple Sources 

IRS 43 (also known as YLW 15) was found to be a binary VLA source with 0'.'6 separation 
(Girart, Rodriguez, & Curiel 2000). IRS 43 is also part of a wide-binary system with GY 263 (see 
Allen et al. 2002 and Table 3). We find IRS 43 to be multiple at 10 /xm (0^'51 separation; P. A. 
= 332.7 degrees) but single in the near-IR. The 2MASS catalog position (detected at H and K) 
for this object is a = 16''27"'26.^94, 5 = -24°40'50" (J2000). This 2MASS position is accurate to 
better than 0'.'2, and it is located very near the VLA 2 component in the 3.6cm map of this system 
presented by Girart et al. (2000). Therefore all the near-IR measurements of this object are likely 
associated only with the south-east AT-band component (see Figure 1). We have imaged IRS 43 
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at 2 /xm wavelength with the 10-m Keck II telescope under conditions of (X'5 seeing but still only 
detected a single source. We estimate that the companion is likely to be i^T > 12 mag from this 
observation and similar ones of others (L. Prato 2002, private communication). 

The brighter mid-IR source in IRS 43 corresponds to VLA 2 and the heavily- veiled, Class 
I near-IR source. VLA 1, the resolved thermal jet component, is at the position of the fainter 
mid-IR source to the northwest. Both VLA/mid-IR sources associated with IRS 43 are embedded 
in extended, faint near-IR nebulosity imaged with HST/NICMOS (Allen et al. 2002). The total 
1.3mm flux from IRS 43 is 75 mJy in an 11" beam (Andre & Montmerle 1994). This level of 
millimeter dust continuum emission is typical of Class I objects in p Ophiuchus, and is therefore 
most likely associated with the heavily-veiled NIR/bright mid-IR source positionally coincident 
with VLA 2. Any protostar at an even earlier evolutionary stage than IRS 43/ VLA 2 would have 
to have an order of magnitude greater millimeter flux than is observed (Andre, Ward-Thompson, 
k. Barsony 1993; Andre, Ward-Thompson, k, Barsony 2000). Given the strength of its mid- infrared 
and sub-millimeter emission, if an embedded, late-stage protostar were present in VLA 1, we 
would have detected it in the near-infrared images as well. The fact that we did not, argues for an 
outflow origin for the mid-infrared emission associated with VLA 1. 

IRS 43 is, in fact, associated with a compact molecular outflow which is nearly pole-on 
(Bontemps et al. 1996). The blue lobe is to the north-west, and the red lobe is displaced south-cast. 
IRS 43 is also a hard X-ray source with strong flares (Grosso et al. 1997), and its total near-to 
far-IR bolometric luminosity is approximately 10 L© (Wilking, Lada, &: Young 1989). Greene 
k, Lada (2000, 2002) did not detect any photospheric absorption features in the high-resolution 
near-IR spectrum of this object, indicating a very large near-IR veiling which is consistent with a 
high mass accretion rate (M = 1.7 x 10~^ yr~^). 

IRS 51 appears to be another example of a flat-spectrum source, in which the extended, 
possibly binary morphologies observed in the near-infrared and in the mid-infrared, do not agree 
with each other (see Figure 2). Although IRS 51 is resolved along an approximately north-south 
direction at both wavelengths, the secondary component at A^-band is at a different position 
angle (PA = 15°) and separation (~ 0'.'7) from the primary than the low-level extension observed 
towards the North in the if -band contours (PA = 10°, separation ~ 1''5). This behaviour suggests, 
analogously to the case of IRS 43, that the mid-infrared "secondary" may very well be a bright 
knot in the molecular outflow from IRS 51 (Bontemps et al. 1996), along a North-South cavity 
seen in scattered light at if-band. 

The near-to-far IR bolometric luminosity of IRS 51 is measured to be only 1.4 Lq. If this 

luminosity is entirely due to a single pre-main-sequence star without accretion, then PMS models 
indicate that its mass would be approximately 0.5 Mq. Greene & Lada (2000) detected weak and 
broad 2.2935 /xm CO band head absorption in the brighter if-band component. This is indicative 
of significant i^C-band veiling (r^ > 1) and fast rotation {v sin z ~ 50 km s~^). Luhman k Rieke 
(1999) determined that the spectral type of this object is K5 - K7, consistent with the low mass 
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suggested from the total system bolometric luminosity. 

The WL 1 system has the SED of a flat spectrum object. It was first resolved as a binary 
system in the near-IR by Ressler (1992), and reported as a mid-IR binary by Barsony & Ressler 
(2000), who classified WL IN as a Class II object and WL IS as a flat spectrum object, based on 
the individual components' spatially-resolved 2-10 micron spectral indices. The WL 1 system has 
a bolometric luminosity of 1.8 Lq (Wilking, Lada, &: Young 1989). It has a peak 1.3 mm flux of 
25 mJy but is not associated with a molecular outflow (Andre Sz Montmerle 1994). Luhman Sz. 
Rieke (1999) determined that the spectral type (of the brighter near-IR component) is later than 
M3, consistent with the low bolometric luminosity for the system. 

The L1689 SN02 system is located in the L1689 dark cloud, southeast of the main L1688 p 
Oph cloud and its embedded stellar cluster. This system is identified as L1689 #5 by Greene et 
al. (1994) who determine that it has a flat 2-10 /^m spectral index (a = 0.0) and a 1.25-20 /um 
calorimetric luminosity lower limit of 1.3 Lq. Bontemps et al. (2001) estimate that its total stellar 
luminosity is 1.9 Lq and its disk luminosity is 1.4 Lq. This source is approximately 2.5' west of 
the peak DCO-I- core position in the southern region of the L1689 cloud. 

Eiroa et al. (1987) flrst found that the Serpens SVS 20 system was binary via CCD imaging 
at 9140A and had extended emission in CCD images. These authors reported a I'.'l separation at 
P. A. = 10°. The binary was subsequently imaged by Huard, Weintraub, & Kastner (1997). These 
authors found flux ratios of SVS 20(S)/SVS 20(N) = 2.89 +/-0.05, 2.99 +/- 0.03, and 2.87 +/- 
0.03 at 3.3, 3.4, and 4.0 microns, respectively. For comparison, our corresponding 10.8 micron flux 
ratio is 2.85. Casali & Eiroa (1996) found that its 2 ^m 5v = 2 CO bands were in emission weakly. 
The 2-14 iim (K band and ISO) spectral index of SVS 20 is found to be flat (a = 0.0) by Kaas 
(2001; private communication). Hurt &: Barsony (1996) found that the far-IR energy distribution 
of this object is also relatively flat over the 12-100 /xm IRAS bands. They estimate that this 
system has 0.3 Mq of circumstellar material. Casali et al. (1993) note that this source has weak 
sub-mm emission, detected at 800 and 1100 /um. They estimate that the bolometric luminosity of 
this system is 0.6 - 30 Lq. By integrating under the spatially resolved SEDs of the SVS 20 system 
presented here for the first time (see Table 3 &; Figure 10), we find L = 4.1 Lq for SVS 20N and L 
= 12.5 Lq for SVS 20S. These luminosity determinations represent only lower limits for the true 
bolometric luminosities of these sources, however, since the integrations were carried out over the 
limited wavelength range (1.25 - 10.8 /xm) over which we have spatially resolved data. 

EC 92 and EC 95 in Serpens form a binary system at 5''03 separation at P. A. ~ 352°. EC 95 
is by far the brighter at 10.8 micron, and is an extremely X-ray luminous proto-Herbig Ae/Be star 
(Preibisch 1999). 

GY 21 and GY 23 (sep'n = 10^'5; P.A. ~ 323°) form a candidate wide binary system. GY 
23 (also known as Elias 23) has been found to be variable by Elias (1978). GY 21 is associated 
with the centimeter source LFAM 3 (Leous et al. 1991)), is classified as a Class II MO source by 
Wilking et al. (2001). 
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Finally, EC 82 (also known as SVS 2) forms a possible wide binary with EC 86 (sepn. = 8''65; 

P. A. 63°), and displays extended nebulosity at J, H, and i^-band. EC 82 is the illuminating 
source of the famous Serpens Reflection Nebula, and is especially apparent in our J and i?-band 
images (see Figure 8). 

4.2. Multiplicity and SED Characteristics 

We define the companion star fraction (CSF) to be the ratio of the number of detected 
binary and multiple systems, to the observed number of single, binary, and multiple systems. For 
reasons discussed in the previous section, we consider IRS 51 to be a single source plus reflection 
nebulosity associated with the blue-shifted outflow lobe from the primary rather than a true 
binary star. Thus, we find a CSF for the Classl/flat-spectrum sources surveyed of 7/19 (37% ib 
14%). This is consistent with the CSFs derived for PMS stars in each of the Taurus, p Ophiuchus, 
Chamaeleon, Lupus, and Corona Australis star-forming regions over a separation range of ~ 150 - 
1800 AU (Leinert et al. 1993, Simon et al. 1995, Ghez et al. 1997, Allen et al. 2002). In contrast, 
the CSF for solar-type main-sequence stars in the solar neighborhood in this separation range 
(11% lb 3%; Duquennoy Mayor 1991, Fischer & Marcy 1992) is approximately one-third that 
of our sample. As is thought to be the case with TTSs, it is possible that this difference is due to 
evolutionary effects or environmental differences. For example, p Ophiuchus and Serpens are dark 
cloud complexes, however some stars form in giant molecular cloud environments (e.g. Lada 1992, 
Carpenter, Snell Sz Schloerb 1995, Phelps Sz Lada 1997). Indeed, the observed binary fraction in 
Orion is consistent with that for solar-type main sequence stars (Prosser et al. 1994, Petr et al. 
1998, Duchene 1999 ). 

Finally, the composite SEDs of SVS 20 and WL 1 are indicative of Class I and fiat spectrum 
sources, respectively. The individual components of SVS 20 both show flat spectrum SEDs, while 
both components of WL 1 appear to be Class II objects. Such behavior is consistent with what 
one typically finds for TTSs, where the companion of a classical TTS also tends to be a classical 
TTS (Prato & Simon 1997). 

5. Summary and Conclusions 

We have obtained new near- and mid-IR observations of 19 Class I/flat-spectrum objects 
in the p Ophiuchi and Serpens dark clouds. These observations are part of a larger systematic 
infrared multiplicity survey of Class I/fiat-spectrum objects in the nearest dark clouds. The 
primary results and conclusions derived from the present study can be summarized as follows: 

1. The CSF for the Class I/flat-spectrum sources surveyed is 7/19 (37% ± 14%). These 
results are consistent with the CSFs derived for PMS stars in each of the Taurus, p Ophiuchus, 
Chamaeleon, Lupus, and Corona Australis star-forming regions over a similar separation range. 
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However, the CSF for solar-type main-sequence stars in the solar neighborhood in this separation 
range (11% ± 3%) is approximately one-third that of our sample, and may be attributed to 
evolutionary effects or environmental differences. 

2. Among the sources which were found to be binary, and for which SEDs of both sources 
could be constructed, we find that the individual components of SVS 20 and WL 1 exhibit the 
same SED classifications. This behaviour is similar to what one typically finds for TTSs, where 
the companion of a classical TTS also tends to be a classical TTS (Prato Sz Simon 1997). 

3. The strength of the mid-infrared and sub-millimctcr emission observed in IRS 43/VLA 1 
suggests that, if an embedded, late-stage protostar were present in this source, we would have 
detected it in our near- infrared images as well; we did not. In addition, the near- and mid-IR 
"secondary" sources observed in IRS 51 are at different position angles and separations from each 
other. We therefore suggest that the companion sources observed in IRS 43 and IRS 51 may 
not be stellar objects, but rather are associated with molecular outflows from the primary YSO, 
although IRS 43 does form a wide binary with GY 263. 
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Fig. 1. — H, K, L, and A/"-band images of the field containing IRS 43. The wide (~ 7'/0 separation; 
see Table 3) binary with GY 263 is apparent in the H, K, and L-bands. Neither component of 
this wide binary was detected at J-band (see Table 2). North is up, and East is to the left in 
each image. Pixel sizes are Gf.'lAS at HKL and 0'.'138 at A^-band, for a total field of view at HKL 
of ~ 9'.'5 X 9^'5, and 8''7 x 8'.'7 at AT-band. Contour levels for the H-hdjad image are 2.5, 5, 10, 
25, and 50% of the peak value; contour levels for the K-band image are 1, 2.5, 5, 10, 25, and 
50% of the peak value; contour levels in the L-band image are 0.5, 1, 2.5, 10, 25, and 50% of the 
peak value. Contour levels in the A^-band image are at 10, 25, and 50% of the peak value. Our 
MIRLIN observations resolved two distinct components at 0'.'51 separation at a position angle of 
332.7 degrees in IRS 43. The southern component coincides positionally with VLA 2; the fainter, 
northern component coincides with VLA 1, the resolved thermal radio jet (Girart, Rodriguez, &; 
Curiel 2000). This component was not resolved in our near-IR images. 

Fig. 2. — J, H, K, L, and A''-band images of IRS 51. The field of view of each image is 
approximately 9'.'5 x 9^'5. North is up, and East is to the left in each image. Contour levels are 25, 
50, and 75% of the peak value in the J-band image; 0.5, 1, 2.5, 5, 15, and 50% of the peak value in 
the H and K-hand images; 1, 2.5, 5, 15, and 50% of the peak in the L-band image; and 10, 25, and 
50% of the peak value in the A^-band image. Although IRS 51 is resolved along an approximately 
north-south direction at both near- and mid-IR wavelengths, the secondary component at A^-band 
is at a different position angle (PA = 15°) and separation (~ 0'.'7) from the primary than the 
low-level extension observed towards the North in the H, K, and L-band contours (PA = 10°, 
separation ~ 1''5). 

Fig. 3. — J, H, K, L, and A^-band images of the WL 1 binary. The field of view of the J HKL 
images is approximately 4'.'6 on a side with 0'.'148 pixels, while the AT-band image is 4'.'3 on a side 
with 0'.'138 pixels. North is up, and East is to the left in each image. Contour levels in the J-band 
image arc at 25, 50, and 75% of the peak value; contour levels in the K, L-band images are 
at 5, 10, 25, and 50% of the respective peak values, and the A^-band image contour levels arc at 
15, 25, and 50% of the peak value. Note that although both sources are about equal brightness at 
J-band, the southern source gradually increases in brightness towards the longer wavelengths, even 
as the northern source progressively dims towards the longer wavelengths. While the combined 
SED for WL 1 implies a fiat-spectrum object Bontemps et al. 2001, the individual components of 
this source exhibit pronouncedly Class II SEDs. 

Fig. 4. — J, H, K, and L-band images of L1689 SN02. Each image is approximately 9'.'5 on a side. 
North is up, and East is to the left in each image. Contour levels are 5, 10, 25, and 50% of the 

peak value in all panels. Although we do not as yet have resolved 10-micron imaging photometry, 
note the differing behavior of the sources with wavelength: although both are nearly equally bright 
at J-band, the southern source dims much faster towards longer wavelengths than its northern 
neighbor. 

Fig. 5. — J, H, K, and L-band images of the wide binary GY21/GY23. Each image is 
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approximately 19" on a side. North is up, and East is to the left in each image. Contour levels for 
the JH and i^-band images are 0.5, 1, 2.5, 5, 10, 25, and 50% of the peak value. Contour levels in 
the L-band image are at 1, 2.5, 5, 10, 25, and 50% of the peak value. 

Fig. 6.— J, H, K, L, and iV-band images of SVS 20. The pixel scale of the JHKL images is 0^'148 
pixel"^, and 0'.'138 pixel ~^ for the A/'-band image. Each image is therefore approximately 4''7 on 
a side. North is up, and East is to the left in each image. Contour levels in the J-band image are 
at 0.5, 1, 2.5, 5, 12, 20, 50, and 75% of the peak value; contour levels in the H and -ftT-band images 
are at 0.5, 1, 2.5, 5, 15, and 50% of the respective peak values. The ii'-band image was smoothed 
with a 3 X 3 boxcar. Contour levels in the L-band image are at 1, 2.5, 5, 15, and 50% of the peak 
value, and the A'"-band image contour levels are at 1, 2.5, 5, 10, 25, and 50% of the peak value. 
The individual components of SVS20 are separated by 1'.'51, at a position angle of 9.9 degrees. 

Fig. 7. — J, H, K, and L-band images of the wide binary EC92/EC95. Each image is approximately 
6'.'7 X 9'.'5. North is up, and East is to the left in each image. Contour levels are 5, 10, 25, and 
50% of the peak value in all panels. 

Fig. 8. — J, iJ, K, and L-band images of the possible wide binary EC82/EC86. Each image is 
approximately 16''3 on a side. North is up, and East is to the left in each image. Contour levels 
are 0.5, 1, 2.5, 5, 10, 25, and 50% of the peak value in all panels. The famous Serpens Reflection 
Nebula is especially apparent in the J and i?-band images. 

Fig. 9. — Spectral Energy Distribution of the p Ophiuchus protostellar binary WL 1. The SED 
of the northern component (N) has been shifted down by 1.0 in the vertical axis for clarity. Both 
components of this source exhibit Class II SEDs. 

Fig. 10. — Spectral Energy Distribution of the Serpens protostellar binary SVS 20. The SED of 
the northern component (N) has been shifted down by 1.0 in the vertical axis for clarity. Both 
components of this source exhibit flat spectrum SEDs. 



Table 1. Summary of Observations 



Source Region RA (2000) Dec (2000) IRTF JHKL UT Date iV-band UT Date/Location 



IRS 43 


p Oph 


16 


27 


26, 


.90 


-24 40 51.50 


2001 


Jul 


09 


1998 Jun 07/Keck II 


GY 263 


p Oph 


16 


27 


26 


.60 


-24 40 45.90 


2001 


Jul 


09 


1998 Jun 07/Keck II 


IRS 67 


p Oph 


16 


32 


01, 


.00 


-24 56 44.00 


2001 


Jul 


09 


Not Observed 


IRS 51 


p Oph 


16 


27 


39, 


.84 


-24 43 16.10 


2001 


Jul 


09 


1997 Jun 25/Palomar 


YLW 16A 


p Oph 


16 


27 


28, 


.00 


-24 39 34.30 


2001 


Jul 


10 


1997 Jun 26/Palomar 


GY 224 


p Oph 


16 


27 


11, 


.17 


-24 40 46.70 


2001 


Jul 


10 


1998 Jun 07/Keck II 


IRS 63 


p Oph 


16 


31 


35 


.53 


-24 01 28.30 


2001 


Jul 


10 


Not Observed 


WL 1 


p Oph 


16 


27 


04 


.13 


-24 28 30.70 


2001 


Jul 


10 


1998 Jun 07/Keck II 


WL 19 


p Oph 


16 


27 


11, 


.74 


-24 38 32.10 


2001 


Jul 


11 


1996 Apr 23/Palomar 


VSSG 17 


p Oph 


16 


27 


30, 


.18 


-24 27 44.30 


2001 


Jul 


11 


Not Observed 


GY 21 


p Oph 


16 


26 


23 


.54 


-24 24 41.50 


2001 


Jul 


11 


1998 Jul 02/Palomar 


GY 23 


p Oph 


16 


26 


24, 


.00 


-24 24 49.86 


2001 


Jul 


11 


1998 Jul 02/Paloniar 


L1689SN02 


p Oph 


16 


31 


52, 


.13 


-24 56 15.20 


2001 


Jul 


11 


Not Observed 


EC 92 


Serpens 


18 


29 


57, 


.75 


01 12 56.96 


2001 


Jul 


09 


1998 Jun 07/Keck II 


EC 95 


Serpens 


18 


29 


57, 


.80 


01 12 52.00 


2001 


Jul 


09 


1998 Jun 07/Keck II 


ISO 159 


Serpens 


18 


29 


32, 


.00 


01 18 42.00 


2001 


Jul 


09 


Not Observed 


SVS 20 


Serpens 


18 


29 


57, 


.70 


01 14 07.00 


2001 


Jul 


09 


1998 Jun 07/Keck II 


EC 53 


Serpens 


18 


29 


51, 


.20 


01 16 42.00 


2001 


Jul 


10 


Not Observed 


EC 82 


Serpens 


18 


29 


56, 


.80 


01 14 46.00 


2001 


Jul 


10 


Not Observed 
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Table 2. JHKL Magnitudes, 10//m flux, Near-IR Colors, and Spectral Indices for p Oph Sources 



Source 






J 




H 




K 


L 


T? ^ 


(J-H) 


(H-K) 


(K-L) 


0^ 


irio o 




1 Q 


.ZD 


1 9 

IZ. 


.00 


Q 


AA 


« 11 


1 ^A 




6.70 




9 fi7 


U.u 


IRS 43 N 


















0.52 












GY 263 


> 


19, 


.29 


14. 


.19 


12. 


.42 


10.71 


0.03 


> 


5.10 


1.77 


1.71 


-0.2 


IRS 67 


> 


20. 


.26 


14. 


.06 


10. 


.34 


7.80 


— - 




6.20 


3.72 


2.54 


0.7^^ 


ixlo 01 O 






fid 

.ou 


1 1 

11. 




O. 


oy 


D.OO 


u. / o 




6.56 


Z.OO 


l.ol 


-u.o 


IRS 51 N 


> 


19. 


.77 


12. 


.65 


11. 


.18 


9.17 


0.48 


> 


7.12 


1.47 


2.01 


1.0 


YLW 16 A 


> 


19. 


.27 


14. 


.72 


10. 


.48 


7.77 


2.65 


> 


4.55 


4.24 


2.71 


1.6 


GY 224 




20. 


.04 


13. 


.33 


10. 


,84 


8.51 


0.25 




6.71 


2.49 


2.33 


0.3 


IRS 63 




16. 


.73 


11. 


.76 


9. 


,34 


7.18 






4.97 


2.42 


2.16 




WL 1 S 




17. 


.21 


12. 


.84 


10. 


,76 


9.02 


0.08 




4.37 


2.08 


1.74 


-0.5 


WL 1 N 




17. 


.23 


12. 


.91 


10. 


,85 


9.18 


0.05 




4.32 


2.06 


1.67 


-0.7 


WL 19 


> 


19. 


.28 


15. 


.02 


11. 


,18 


8.32 


0.18 


> 


4.26 


3.84 


2.86 


0.2 


VSSG 17 




17. 


.05 


11. 


.34 


9. 


,07 


6.85 






5.71 


2.27 


2.22 


0.3^= 


GY 21 




14. 


.59 


11. 


.87 


10. 


,09 


8.04 


0.40 




2.72 


1.78 


2.05 


0.1 


GY 23 




10. 


.74 


8, 


.59 


7. 


,39 


5.92 


2.02 




2.15 


1.20 


1.47 


-0.4 


L1689 SN02 S 




12. 


.13 


10. 


.05 


8. 


,86 


7.96 






2.08 


1.19 


0.90 




L1689 SN02 N 




12. 


.16 


9 


.81 


8. 


,32 


6.64 






2.35 


1.49 


1.68 





''lO /um fluxes are listed in Janskys. 

''Spectral index not calculated if both MIRLIN and ISO (Bontemps et al. 2001) fluxes 
unavailable. 

"^Spectral index calculated using the 6.7 /im and 14.3 \jxa. ISO fluxes from Bontemps et al. 2001. 
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Table 3. JHKL Magnitudes, lO/xm flux, Near-IR Colors, and Spectral Indices for Serpens 

Sources 



Source 


J 


H 


K 


L 




(J-H) 


(H-K) 


(K-L) 


a*- 


EC 92 
EC 95 


14.71 
15.93 


12.03 
11.90 


10.44 
9.78 


8.96 
8.34 


0.55 
0.10 


2.68 
4.03 


1.59 
2.12 


1.48 
1.44 


0.7 
-0.4 


ISO 159 


11.25 


9.30 


8.33 


6.68 




1.95 


0.97 


1.65 




SVS 20 S 
SVS 20 N 


12.20 
13.96 


8.74 
10.30 


7.09 
8.38 


5.34 
6.29 


4.36 
1.53 


3.46 
3.66 


1.65 
1.92 


1.75 
2.09 


-0.1 
0.0 


EC 53 


> 19.31 


16.32 


14.65 


11.30 




> 2.99 


1.67 


3.35 




EC 82 
EC 86 


11.75 
12.85 


10.32 
11.57 


8.99 
10.90 


7.24 
10.49 




1.43 
1.95 


1.33 
0.67 


1.75 
0.41 





'^lO fluxes are listed in Janskys. 

''No spectral index calculated if 10 /xm flux unavailable. 
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Table 4. Separations and Position Angles for the Binary Sources 



Source 


Separation (") 


Position Angle (degrees)'' 


IRS 43/GY 263 


6.99 


322.0 


WL 1 


0.82 


321.2 


GY 23/GY 21 


10.47 


322.6 


L1689 SN02 


2.92 


240.3 


SVS 20 


1.51 


9.9 


EC 95/EC 92 


5.03 


352.1 



''Measured with respect to the brightest source at iiT-band; 
the southernmost source in all cases. 
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IRS 43 





GY263 



K 






N 



IRS 43 
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IRS 51 

T 
J 
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T I 

ri 


p 
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IN 
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WLl 

T 
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T I 

ri 


W) 




❖ 
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IN 
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L1689 SN02 
J 




H 




K 






L 
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1 . • 

J 

- •'-^ (. 
' \ / 


^^^^^^ * - " 
♦ 


GY21 
Elias 23 (GY23) 


K 
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J 

EC92 
EC95 


H 


K 

^) 
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2 1 1 1 1 1 1 1 1 1 1 r 

: SVS20 




0.5 

log A (yu.m) 



